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ABSTRACT 



Shear flow instabilities can profoundly affect the diffusion of momentum in 
jets, stars, and disks. The Richardson criterion gives a sufficient condition for in- 
stability of a shear flow in a stratified medium. The velocity gradient V' can only 
destabilize a stably stratified medium with squared Brunt- Vaisala frequency A^^ 
if V'"^ /4: > N"^. We find this is no longer true when the medium is a magnetized 
plasma. We investigate the effect of stable stratification on magnetic field and 
velocity profiles unstable to magneto-shear instabilities, i.e., instabilities which 
require the presence of both ma gnetic field and shear flow. We show that a family 



of profiles originally studied by 



Tatsuno fc Dorland 



( 120061 ) remain unstable even 



when < A^^, violating the Richardson criterion. However, not all magnetic 

fields can result in a violation of the Richardson criterion. We consider a class 



of flows originally considered by 



Kent 



( 119681 ). which are destabilized by a con- 
stant magnetic fleld, and show that they become stable when V"'^/4 < A^^, as 
predicted by the Richardson criterion. This suggests that magnetic free energy is 
required to violate the Richardson criterion. This work implies that the Richard- 
son criterion cannot be used when evaluating the ideal stability of a sheared, 
stably stratifled, and magnetized plasma. We briefly discuss the implications for 
astrophysical systems. 



Subject headings: instabilities; MHD; stars: magnetic fields; stars: rotation; Sun: 
magnetic fields; Sun: rotation 
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1. Introduction 

Rotation plays an important role in the structure and evolution of stars. Although 
rotation directly modifies hydrostatic equilibrium only in the most rapid rotators, it drives 
large scale circulation, modifies the structure of convection and the nature of convective 
transport, and is a key component of magnetic dynamos. These phenomena in turn modify 
the rotation through a complex interplay of nonlinear processes. 

Shear fiow instability is one of the mechanisms through which rotation influences and 
is influenced by its environment. The motion associated with the instability generates 
stresses, which react back on the flow and drive it toward a stable state. If the amplitude of 
the unstable perturbations is sufficiently large, the motions become turbulent. Shear fiow 
instability and shear fiow turbulence can amplify magnetic fields and mix chemical species, 
in addition to modifying the rotation profile itself. 

In the case of the Sun, and possibly other low mass main sequence stars, the most 
likely venue for shear fiow instability is the so-c alled tachoclin e, the region of strong shear 
just below the base of the convection zone (see 



Gough 



20071 . for a review). Although 



the mechanisms which maintain the tachocline are still uncertain, it is almost certainly 
a component of the solar dynamo, and its existence has implications for the way the 
convection zone, which is spun down by the solar wind, is coupled t o the radiativ e core. 



The tachocline may be 



of the field ( 



( lOgilvie 



subject to purely hydrodynamic instabilities ( jRashid et al. 



Oilman fc Foxl ( 119971 ) 



instabilities ([Silvers et al. 



mm, 



Kitchatinov &: Riidigerl (120091)'). gl o bal MHD instabilitie s driven by the latitudinal structure 



Gilman et al 



( 120071 )) magnetorotational instabilities 



2007[). and, if hydroma g netic forces are la r ge en ough, magnetic buoyancy 



(120091 ). IVasil fc Brummelll t00% ). All these instabilities could 



modify the tachocline's structure. 



Massive stars, which evolve quickly and tend to rotate rapidly, are potentially more 
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profoundly affected by shear ffow instability. The past two decades have witnessed 
significant advances in understanding how the inter nal rotation of massive lum inous 
stars shapes, and is shaped by, their evolution (see 



Maeder &: Meynet 



2000a. and 



references therein for a comprehensive review). Rapidly rotating massive stars follow 
bluer, more-luminous evolutionary tracks in the Hertzsprung-Russell diagram (HRD) than 
non-rotating equivalents, because strong meridional circula tion injects fresh hydrogen fuel 



into the convective core (see, e.g. 



Meynet fc Maeder 



20001 ). This rotational mixing brings 



CNO-cycle nucleosynthetic products from the stars' cores to the ir surfaces, leading to 
changes in photospheric abundance ratios (e.g.. lTalon et al.lll997l ). 



The prev ailing view o f rotation in massive stars is based on a canonical narrative 



developed by IZahnI ( 1l992l ). In this scenario, turbulent diffusion of angular momentum 
is highly anisotropic, with much stronger transport in the horizontal direction than the 
radial one. This leads to a 'shellular' rotation profile, in which the angular velocity is 
constant on spherical shells. The exchange of angular momentum between these shells is 
then mediated by a combination of meridional circulation, convection (in convective zones) 



and radial turbulent diffus ion. The turbulence itself is driven by secular s^ 



(Maeder & Mevnet 


2000c 


i), which grow 


s on 


Maeder & Mevnet 


1996 




Talon & Zahn 


1997) 



lear instabilit y 



Maeder 



1995 



Recent studies have considered the role that magnetic f ields might play in modifying 



angular momentum transport (e.g. 



Maeder fc Meynet 



20041 ). Generally, these studies of 



the impact of magnetic fields have focused around co ntributions to the ra dial angular 



m omentum diffu sivity arising from the field stiffness (jPetrovic et al. 



as 



Spruit 



20051 ). However, 



(119991 ) has discussed, a field can also introduce new instabilities that play a 
role in angular momentum transport. In this paper we explore a hitherto-overlooked 
magnetic-mediated instability, whereby the presence of a horizontal field can destabilize a 
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stratified shear layer that — according to the Richardson criterion — would otherwise be 
stable. 

The paper is organized as follows. First we will briefly discuss shear flow instabilities 
in ^ In ^ we set up the eigenvalue problem which determines the linear stability of 
an MHD shear flow in a stratified medium. We review previous analytic results in ^ 
and describe our numerical methods for solving the eigenvalue problem in ^ Starting 
in ^ we examine specific examples, first adding stratification to the linear velocity and 
parabolic magnetic field example considered in a recent paper by Tatsuno & Borland (2006) 
(hereafter TD06). Our key result is that sufficiently strong parabolic magnetic fields can 
yield instability for arbitrarily strong stratification, in violation of the Richardson criterion. 
We consider and extend a family of velocity profiles which Kent (1968) (hereafter K68) 
showed can be destabilized by a constant magnetic field in ^ Contrary to the parabolic 
magnetic field case, it seems that the introduction of a constant magnetic field cannot 
result in a violation of the Richardson criterion. This suggests that the free energy of an 
inhomogeneous magnetic field is essential to breaking the Richardson criterion. We discuss 
possible applications to rotating stars in ^ and conclude in ^ 



Introduction to Shear Flow Instabilities 



The best known shear flow instability is the hydrodynamic Kelvin- Helmholtz instability. 
The Kelvin-Helmholtz instability has been studied extensively. Perhaps the most famous 
result is the inflexion point criterion, stating that a necessary condit ion for instabili t y is t he 



presence of an inflexion point in the velocity profile (see, for example. 



Drazin fc Reidl fllQSlI )). 



Others have als o given necessary condit i ons for instability, making extra assumptions on 



the flow profile (ILin 



1955 



Howard 



1961 



Rosenbluth fc Simon 



19M). 
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Many have worked to extend parts of these results to magnetohydrodynamic shear 
instabihties. It is well know n that a sufficiently s trong magnetic field stabilizes the 



Kelvin-Helmholtz instability (j Chandr asekhar 



196l[ ). It was shown years ago, but is perhaps 



less well known, that a magnetic field can destabilize an otherwise stable shear fiow 
(K68). In particular, an infiexion point is no longer necessary for shear instability. In 
hydrodynamics vorticity is froze n into the fiow, ensuring that perturbations are stable 



when there is no infiexion point ((Lin 



19551 ). but the presence of a magnetic field can break 
the vorticity frozen-in condition, relaxing the infiexion point criterion. K68 constructed a 
family of fiow profiles which are marginally stable in the absence of a magnetic field and 
destabilized by a uniform field parallel to the direction of fiow. TD06 studied how a linear 
fiow profile, which has no infiexion point and is marginally stable, can be destabilized by 
a particular family of magnetic field profiles. In particular, TD06 find that a parabolic 
magnetic field can render a linear velocity profile unstable. 

In this paper, we add a new piece of physics to the analysis: density stratification. We 
employ the Boussinesq approximation and assume that the plasma is stably stratified, i.e. 
the squared Brunt- Vaisala frequency, A^^, is positive. In hydrodynamics, the Richardson 
criterion provides a suffici ent condition fo r the stability of a shear fiow in a stratified 



medium (see, for example. 



Drazin fc ReidI ( Il98ll )). The interchange of two fiuid elements 



at different heights can release kinetic energy from the flow. A necessary condition for 
instability is that the gravitational energy required for the interchange must be less than 
the kinetic energy released. However, in the presence of an inhomogeneous magnetic fleld, 
energy can also be extracted from the magnetic fleld, even if the fleld would be stable in the 
absence of shear flow. Our main result is that the Richardson criterion no longer holds for 
inhomogeneous magnetic flelds. 



We will only consider the effect of stable stratiflcation on magneto-shear instabilities. 
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However, iTatsuno et al.l (|2003[ ) studied how a shear flow can destabihze a homogeneous 
magnetic field in the presence of an unstable density gradient. They found that a linear 
(Couette) velocity profile can be destabilizing when the velocity shear was not too strong. 
Their result is similar to our own in the sense than the system is maximally destabilized 
when the velocity gradient, magnetic field, and density stratification all have comparable 
strength. 

In this paper we consider only ideal instabilities, i.e. we set the resistive, viscous, 
and thermal diffusivities to zero. Diffusive effects could unleash a host of additional 
instabilities such as tear ing modes (e.g. 



Furth et al. 



Schmitt fc Rosner 



1963[ ). d oubly diffusive modes (e.g . 



Maeder fc Meynet 



2mm- 



19831 ). and secular shear instabilities (e.g.. 
Although such instabilities are important in their own right, in this paper we focus entirely 
on dynamical instabilities. 



3. Basic Equations 

The time evolution of an ideal, incompressible plasma is given by 

^ + V.Vvj = -v(p+-)+-B.VB-9Pe.. (1) 

V X (V X B) , (2) 



dB 

'dt 

= V-V, (3) 



= V-B, (4) 
= f + V-Vp, (5) 

where the symbols have their usual meanings. Equation ([T]) is the momentum equation, 
eqn. ([2]) is the induction equation, eqn. ^ enforces incompressibility, eqn. (jl]) is the 
divergenceless magnetic field condition, and eqn. ([5]) is the continuity equation. We will 
write the unit vectors in the x, and z directions as e^, e^, and respectively. The 
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gravitational strength is parameterized by g, and gravity is assumed to point in the — 
direction. We denote background velocity and magnetic fields with capital letters, and 
then perturb the background fields with fields denoted with lower case letters, except that 
the background density is denoted p, and the perturbed density p. We assume that the 
background quantities p, V, B are all functions of only z, and that our domain is the volume 
between z = —zq and z = +zq with "free-slip," perfectly-conducting boundary conditions in 
the z direction, and periodic boundary conditions in the x and y directions. By "free-slip," 
we mean no constraint on perturbed quantities in the x and y directions at the walls, but 
that perturbations have no z component at the walls. These are the boundary conditions 
adopted by TD06 (who termed them "no-slip" which is not correct — as will be shown in 
§6.3[ the perturbations slip along, but do not penetrate, the walls). Next, we assume that 
V is oriented in only one direction throughout the domain, which we define to be the x 
direction. Thus, we take 



in Cartesian coordinates. The background fields are assumed to be in equilibrium, so we 
have that 



Equation ([8]) specifies an integral equation for the background pressure p for arbitrary 
B and p. The induction equation and continuity equation for the background fields are 
automatically satisfied by the geometry we have imposed. 

Now assume the perturbation fields all have the form 



V=(\/(z),0,0), 



(6) 



in Cartesian coordinates. The background magnetic field B is 



B = {B,{z),By{z),0) 



(7) 





/(x, y, z, t) = f{z) exp{ik^x + ikyy - ik^ct). 



(9) 
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We will take k = k^e^ + kySy, k = |k| and k = k/fc. In many applications, the density 
gradient p' is small in comparison to the velocity gradient V — where prime denotes 
differentiation with respect to z — but the strength of gravity g is large. Assuming this, we 
recover the Boussinesq approximation, in which we drop terms proportional to p' alone, but 
keep terms proportional to gp'. These assumptions yield the following eigenvalue problem 
for ^, the plasma displacement in the z direction: 

( [kl {V - cf - k'A^] e)' - e [kl {V - cf - k'A'] i + k'N^i = 0, (10) 

where A = k • 'B/^/ppo is the Alfven velocity, and A^^ = gp'/p is the Brunt- Vaisala 
frequency in the Boussinesq approximation. To simplify our analysis, we assume that A^^ is 
constant throughout the domain, which corresponds to the exponentially decaying density 
profile. When computing the Alfven velocity, the Boussinesq approximation will allow us 
to consider p to be a constant. The boundary conditions are that = at the boundaries 
at z = —zq and z = +zq. 

There is an asymmetry in how velocity shear, magnetic fields, and density stratification 
depend on the wavenumber k. For k = kyBy, k^ = and the velocity shear is irrelevant 
(note that kxC, the growth rate, could still be finite). The purpose of this paper is to 
examine the interplay between velocity and magnetic fields, so we will not consider this 
case. Also note that the Alfven velocity, as it occurs in eqn. flTOj) . is a function of k. 
For example, if B is constant in the z direction, there exists a k for which A = 0, so the 
magnetic field would have no effect on such a perturbation. The strength of gravity in 
relation to shear flow contains a factor of k'^/k'^. Thus, gravity is maximally destabilized by 
shear flows when ky = 0. 

Consider an eigenvalue problem for magnetic field B, velocity V, Brunt- Vaisala 
frequency A^^, and wavenumber k = fc^^e^. + kyCy, with ky ^ 0. We will show that this 
eigenvalue problem is equivalent to another eigenvalue problem with ky = 0, but with 
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different B, A^^, and K. Define B' = e^k ■ B/k^, N''^ = ^N^/kl and k' = ke^. Then 
the magnetic field B', velocity V, Brunt- Vaisala frequency A^'^, and wavenumber k' 
have the same eigenvalue equation as above. Thus, finite ky is equivalent to ky = 0, if 
one appropriately rotates and augments the magnetic field, and increases the density 
stratification. With this in mind, we will consider the ky = case in the remainder of this 
paper, for which the eigenvalue equation reduces to 



( [{V - cf - A'] 0' - e [{V - cf -A^]i + N'i = 0. (11) 

The eigenvalue eqn. ( ITT]) possesses some symmetries. First, the sign of A is 
unimportant, so changing the sign of the magnetic field does not change the problem. 
Another symmetry is translational: taking V ^ V + AV and c — ?■ c — AV corresponds to 
Galilean transformations. Thus, without loss of generality, we can and do put ourselves 
in a frame in which V{0) = 0. To make the problem more tractable, we add additional 
symmetries to the equation by postulating that A is even in z and V is odd. There is a 
rescaling symmetry: eqn. (ITTi) remains invariant under 

Z z/zq, 

V ^ V/zo, 

A ^ A/zo, (12) 
k — )■ kzQ, 
c — )■ c/zq. 

Note that A^^ and Ri are left unchanged under this transformation. 

There is also structure in the eigenvalues. In general, c and ^ are complex; c = Cr + ici, 
^ = + i^i- We ignore the singular Cj = case. If is an eigenfunction with eigenvalue c, 
then the complex conjugate of ^c, is a solution to eqn. ( fTTl) with eigenvalue c*. Thus, 
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eigenvalues come in complex conjugate pairs, regardless of the symmetry properties of A 
and V. Assuming that A is even and V is odd, we can show that if c is an eigenvalue, then 
— c is also an eigenvalue, with eigenfunction C,c{—z)- 

Numerically, we only find eigenvalues with Cr = 0, and with the following eigenfunction 
symmetry: If we normalize the eigenfunction ^ such that ^(0) = 1, then C,r is even and C,i is 
odd. In §^ and Ul we assume that = and the eigenfunction has this symmetry. These 
properties are linked. If we multiply eqn. ( 1111) by ^* and integrate over the domain the 
result is 



[iv - cr - A^] ( 

J -20 y 



dz 



+ \-N'\i\'dz = Q. (13) 



The imaginary part of eqn. (1131) is 



2lCi / {Cr-V) 



20 



di ' 



dz 



+ k'^\i\^\dz = Q. (14) 



If the real and imaginary parts of each have definite parity, the term proportional to V 
in eqn. f lT^ vanishes. Therefore Cr-Cj = 0, and unstable modes have Cr = 0. This result is 
useful in searching for unstable modes, as described in ^ 

We find that generally the growth rate c = is small in comparison to V, which 
is C(l). When V'^ = A^, the coefficient of the ^" term in eqn. ( ITTj) goes to | — c^^l <^ 1. 
Thus, the equation becomes "almost singular" when \V\ = \A\, and becomes actually 
singular when c = 0. The "almost singularities" are characterized by large gradients in the 
eigenfunctions, as is shown in §^ and [3 

We will often consider the limit k"^ = 0. When k"^ = 0, the growth rate, kc, is formally 
zero. However, one can view the eigenvalue c as a function of the various parameters 
A, V, k"^, N"^. We assume that c(A;^) is analytic about k"^ = 0, so our results for the k"^ = 
case still hold in a neighborhood of k"^ = 0. Thus, when we consider k"^ = 0, we are really 
taking the limit as k becomes small. The k"^ term in eqn. ffTTj) is only important when it is 



comparable to the scale heights of the velocity and magnetic fields and the perturbation ^. 
Numerically, we find that kzo < 0.1 is "small" for the examples presented in this paper. 



4. Review of Analytic Results 

Shear fiow instabilities are global instabilities. Thus, the two categories of analytic 
results — necessary conditions for instability and sufficient conditions for instability — can 
be viewed as local and global conditions. Necessary conditions for instability give criteria 
which must be satisfied in at least one spot in the domain, whereas the sufficient conditions 
for instability are global criteria involving integrals over the domain. We present a short 
overview of the analytic results regarding the linear stability of shear fiows. We begin by 
discussing shear fiows alone, and then add stratification, a magnetic field, and then both. 
The zero magnetic field and zero density gradient cases can be viewed as limits of the more 
general problem. 



4.1. Shear Flow Instabilities 



Probably the best known result is the infiexion point criterion, which states that 
V" must have a zero in the domain for there to be instability. This is a local, necessary 
condition. There are several physical interpretations of the infiexion point criterion. 
Consider the Reynolds stress of the perturbation, r = —pv^Vz, where the bar denotes 
averaging with respect to x. A s sumin g c 7^ 0, one can show that dr/dz has a zero iff V" 
has a zero (for instance, in iLinI (119551) o r K68). Since r = at the boundaries, when c 7^ 0, 



we must have that V" has a zero. 



Linl (119551 ) has proposed an alternate interpretation 



considering vorticity. A zero in V" corresponds to an extremum in vorticity, and Lin has 
shown that perturbations feel a restoring force unless they are at an extremum of vorticity. 
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The inflexion point theorem is useful because it rules out a large class of velocity 
profiles as stable. However, it c annot be used to show that a particular shear flow is 
unstable. 



Rosenbluth &: Simod (119641 ) were able to prove a necessary and sufficient 
condition for instability by using the additional assumptions that V" has a single zero and 
V is monotonic. Under these assumptions, V is unstable in zi < z < Z2 if and only if 



V'{V, - V) 



22 



V" 



-dz > 0, 



(15) 



v'^iv-v,: 

where Vc is the velocity at the inflexion point. This result is derived for the fc^ = case. A 
priori, it seems that there could be velocity profiles which are unstable for k'^ > but stable 
for k"^ = 0. Then an instability co nditi o n for Z;;^ = would be only sufficient for instability. 



This is addressed by a theorem of iLinI (jl955l ) which shows that under the assumptions of 
Rosenbluth & Simon, velocity profiles which are unstable for k^ > are also unstable for 
P = 0. 



4.2. Shear Flow Instabilities in a Stratified Medium 



The key stability result for stratified media is the Richardson criterion, a necessary 
condition for the instability of a shear flow in a stratified medium. If 



1 



(16) 



everywhere, then the re i s stability. A 



Chandrasekharl (jl96ll ) or 



p hysica l interpretation (see, for example. 



Drazin &: ReidI (Il98ll )) is that if exchanging fluid elements at 



slightly different heights increases the potential energy more than it decreases the kinetic 
energy, then the perturbation is stable. 



Provided that Ri < 1/4, we have that 



fc^c? < max ( ^V'^ - 



(17) 
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discussed in 



Drazin fc Reidl f|l98lh ) 



This result bvlHowardI (11961) follo ws from the proof of the Richardson criterion and is also 



4.3. Magneto-Shear Instabilities 

Magnetic fields can both stabilize and destabilize shear flows. First we consider their 
stabilizing effect. Perturbations which bend magnetic fieldlines induce a restoring magnetic 
tension force. A classic result is that in a constant density medium, the vortex sheet 
V{z) = —U for z < Q and V{z) = +U f or z > for some con stant U, is stabilized by a 



magnetic field A if and only if > V"^ (j Chandr asekhar 



196l|). This 



step function velocity 



Keppens et al 



mm 



profile is the limiting distribution of V{z) = Uotanh.{z/a) as a — )■ 0. 
have investigated the hyperbolic tangent V case with a constant magnetic field, including 
compressibility, and found the magnetic field stabilizing. These results were qualitatively 
similar to those by Chandr asekhar, which is expected because a constant magnetic field has 
no length scale (or it has an infinite length scale), so it cannot tell the difference between 
the a — )• and a finite case. 

Keppens et al also found that the addition of a non-uniform magnetic field could be 
destabilizing. When they added a small field A{z) = —Aq for z < and A{z) = Aq for 
z > 0, they found that the growth rate increased, and was even larger when A reversed 
smoothly. Although their calculation, unlike ours, includes compressibility, there is one 
robust effect which is always present: magnetic fields allow transfer of vorticity between 
fiuid elements. The loss of the frozen-in vorticity constraint changes the range of motions 
allowed in the plasma, and yielding instability. 

We now review some general results on magnetoshear instabilities in order to 
understand how the Richardson criterion can be violated by the introduction of a magnetic 
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field. 



The necessary and sufficient instability condition of iRosenbluth &: Simon! (Il964r) (eqn 



Chen fc Morrison! ( 



( IT5|) ) has been generalized to the MHD case by K68 and 
arguments use that when k"^ = 0, there is an exact solution to eqn. (fTTjl . 

dz' 



19911). Both 



e(^) 



and then define 



/(c) 



dz 



e(^2). 



(19) 



^ {V-c^-A^ 

The eigenvalues of eqn. fill I) are then just the zeros of /(c), and one can search for 
instabilities by implementing Nyquist's method to determine if there are any zeros of /(c) 
for Cj > . Nyq uist's method is an application of the Argument Principle (see, for instance 



Gamelin! ( 12001! )). which states that the integral of the argument of /(c) on the boundary 



dD of some region D is equal to 27r(A^o ~ A^oo), where A^o is the number of zeros of /(c) 
in D and A^oo is the number of poles of /(c) in D. In our case, we assume A^oo = 0, 
so counting the number of times /(c) wraps around the origin tells us how many zeros, 
i.e. unstable modes , ther e are. Further discussion of Nyquist's method can be found in 



Krall fc Trivelpiecel (11973! ) 



Nyquist's method can only be applied if we know what contour to use. The real part of 
c can b e bounded by extending an important hydrodynamic result by Rayleigh. It can be 



shown (iHughes fc Tobias 



20011 1 that Cr must lie in the range of V, so the contour in c space 



is bounded by Vmin < Cr < Vmax- The lower bound for cy is 0"*", and the upper bound can be 



recovered by modifying Howard's semic ircle theorem ( 



case, Howard showed (see, for instance. 



loward 



Drazin &: Reid 



196l! ). In the hydrodynamic 



(jl98l|)) that 



Cr 2 {.^max ~l~ ^min) 



Thus, we have that c^ < l/2{Vmax - Vmin)- 



+ c < 



n 2 



(20) 



Hughes &: Tobias! (l200l[ ) have shown that in 



- 16 - 



MHD, we have the two inequahties 



and 



1 



+ C < 



1 



\ In 



(21) 



(22) 



This gives an even stronger upper bound on q, that 



Q ^ \/(l/2( V^nax ^min))'^ (^^)mm 



(23) 



These two inequahties can be used to show stabihty, if one can show that there are no c 
which simultaneously satisfy both inequalities. 



Chen fc Morrison! (jl99ll ) used Nyquist's method to provide a sufficient condition for 
instability for flows in which V is even, and A is either odd or even. They showed that 

dz 



20 



> 



(24) 



as e — )■ is sufficient for instability. Note that it is not assumed that V has an inflexion 
point. 

K68 considered the effects of a small, constant magnetic field on a stable velocity 
profile. He showed that when V" has a single zero, and there exist points ysiUt such that 
the velocities at these points, Vs, Vt satisfy Vs — Vt = 2A and — VI = 0, then 

dz 



M{A) ^ p 



> 0, 



(25) 



implies instability. Here, Cq is defined by Cq = (Vg + Vt)/2, and p denotes the principal value 
of the integral. For small A, cq is the velocity at the infiexion point, but as A increases, it 
can deviate somewhat. For a marginally stable velocity profile, we have M(0) = 0. In the 
remainder of this section, we will use ' (dot) to denote derivative with respect to A. In the 
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limit A — >■ 0, we have M{A) — >■ 0. Thus, to evaluate the stabihty of V to infinitely small A, 
we need to consider M(0), which Kent shows is given by 



M(0) = TcM I , + 2 / , 



where Vq is the velocity at the inflexion point and 



(26) 



c.(0) 



y 



(4) 



(3)- 



(27) 



This criterion is useful because one can change variables to integrate over V ^ and if = 
and ijj{y) := dz/dV is even, then 

"^^ udy 



M(0) 



^4 ' 



(28) 



where = V{zi). Although these conditions are sufficient for instability, they are not 
necessary. Unlike in the hydrodynamic case, there can be unstable modes for finite k"^ for a 
velocity profile which is stable at /c^ = (K68). 

Another way to tackle the general problem with arbitrary velocity and magnetic field 
profiles is to attempt to extend the physical arguments behind the infiexion point criterion 
to the MHD problem. In the MHD problem, one must consider both the Reynolds and 
Maxwell stresses, so the total stress is given by 



(29) 



A necessary condition for instability is still drtot/dz — somewhere in the fiow. K68 has 
shown that this condition can be written as 



1X1 



.X'' 



0, 



or 



'2XX" - X 
AX^ 



/2" 



0, 



(30) 



(31) 
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where X = {V — c)'^ — A? . Unfortunately, these (equivalent) conditions are not as useful as 
the inflexion point criterion because they depend on both the flow profile and the growth 
rate. Thus, one needs to check eqns. fl5Ul) or fl?Il) for all possible c. This condition seems to 
be fairly weak, and is satisfied by many stable profiles. 



4.4. Magneto-Shear Instabilities in a Stratified Medium 

The addition of a magnetic field to a shear flow in a stratified medium makes the 
problem significantly more complex. The Richardson criterion is no longer valid, but it 
can be generalized. We have carried out the same analysis used to derive the Richardson 
criterion, but included magnetic fields. The result is that if 

everywhere in the domain, then the system is stable. Here, Z =\ — A? j{y — cf' . Similarly 
to the generalization of the inflexion point criterion (eqn. (ED]), (EI])), this condition involves 
c. This condition also seems to be weak. 

Although we normally assume t hat Cr = 0, this co n dition can be relaxed, and we can 



flnd bounds for c^. The argument by 



1 r 



Hughes fc TobiasI (l2001|) mentioned in S3] still holds 



when stratiflcation is introduced, and shows that Cr must lie within the range of V . This 
bound on Cr is valid with and without magnetic fleld, and with and without stratiflcation. 



5. Numerical Methods 

Because the problem is global, analytic results exist only in cases with particular 
symmetries (i.e., /c^ = or A^^ = 0), so we must generally solve for stability numerically. 
We have implemented three numerical methods for solving the eigenvalue problem, eqn. 



- 19 - 



f fTTj) . In the first, we discretize the equation onto a Chebyshev grid, and use a finite 
dimensional approximation for the differential operator. Then eqn. f ill I) can be rewritten 
as a generalized finite-dimensional eigenvalue equation: 

( V o \ f v,\ ( -tk^VV + it^V" ikAV - ikA" -N^P \ ( ^ 
7 10 b, = ikA -ik^V b, (33) 

\001/\p/\ 1 -i/c^y I \ P ) 

where V = dl - k"^ . Matlab was used to solve this finite dimensional eigenvalue problem. 
This approach was useful when we did not require high resolution. This method was not 
able to resolve the large gradients in the eigenfunctions that sometimes appeared when 
1^1 = 1^1- 

Another strategy, for = 0, was implementing Nyquist's method. We used 
Mathematica to calculate /(c), as defined in eqn. ( |T9i) for various c. As mentioned in 
§4.41 we know that Cr lies between the minimum and maximum of V . The advantage of 
Nyquist's method is that we need not assume that c is imaginary. We picked the rectangle 
with vertices at ie + Vmax, + Knm, ici + Knm, and ia + Vmax as the contour, with a of order 
one and e small. If one plots /(c) where c traverses this contour, it is easy to see if there are 
any unstable modes with c in this contour. We varied the size of the rectangular contour to 
find the exact eigenvalues. For the examples presented below in §^and[7l eigenvalues were 
always purely imaginary, and the eigenfunctions had the symmetry properties described in 



Finally, we used a finite difference relaxation code to integrate across the domain. We 
assumed that c was imaginary, and integrated eqn. (ITTi) over the domain for c between 
ie and ia for a of order one and e small, in logarithmic steps. When the real part of /(c) 
changed sign between two consecutive steps, the secant method was used to find the zero 
in the real part of /(c), which corresponds to a zero in /(c). This algorithm was the 
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most efficient, but makes the assumption that the eigenvalues are purely imaginary. As 
mentioned in §31 we have not found any eigenvalues with non-vanishing real part using the 
other two methods mentioned above, so this seems to be a valid assumption. 

All three numerical methods give similar results in cases where we used more than one. 



6. Linear V, Parabolic A 

In this section, we add density stratification to the linear velocity and parabolic 
magnetic field profiles considered by TD06. The main result is that we find instability 
even when V'^/A < N"^ everywhere, i.e., when the Richardson criterion predicts stability. 
We believe this is because the magnetic field provides another free energy source for the 
instability. At k"^ = 0, there are magnetic field profiles which are unstable for arbitrarily 
large A^^, but when fc^ > 0, there is only a finite range of A^^ which are unstable for the 
profiles considered here. 



6.1. The Field and Flow Profiles 

Consider the following velocity and magnetic field profiles in a domain from z = —1 to 

z = +l: 

V{z) = z, (34) 
A{z) = {l-a)z'^ + a. (35) 

These are the fields considered by TD06 in Section Ill.A.l (where we call their ai parameter 
a). The magnetic field is a parabola with A{0) = a and A = 1 at the boundaries. 

An important characteristic of these profiles is that neither the magnetic field nor 
the velocity profile are unstable by themselves. The instability is truly a magneto-shear 
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instability, as both magnetic field and shear fiow play a part in rendering the profiles 
unstable. In this respect, this example is different f rom those considere d by others in 



which a magnetic instability is st abilized by grav ity (jPikpati et al. 



destabilizes a stratified medium (Newcomb 



a buoyancy instability ( iHowes et al. 



20011). 



1961 



20091 ) ■ a magnetic layer 



or magnetic field and shear fiow modify 



These profiles can be viewed as local approximations to a wide range of field and fiow 
profiles. The parabolic magnetic field profile is valid locally whenever B has an extremum, 
which we take to be at 2; = 0. As mentioned in ^ taking A — —A does not change the 
problem, so although we are considering a local minimum, the exact same results hold for 
A{z) = — (1 — a)z'^ — a, which characterizes a local maximum. We can always transform 
to a frame in which ^(0) = 0, so the velocity has a local expansion of the form of equation 



To view these profiles as a local approximation, we also need to make an assumption 
about the relative strength and scale of variation of the magnetic field and the shear fiow, 
since we require that \V\ = \A\ at the boundary. When a is close to zero, the magnetic field 
and velocity are changing at similar rates, so the locality assumption is plausible; but when 
a is close to one or very negative, the scale heights of the fiow and magnetic field are very 
different, so viewing these profiles as a local expansion is not as accurate. 

Depending on the sign of a, the magnetic field has either two or zero nulls. When 



a < 0, A = at 



z = ± 



a 



a 



(36) 



When a > 0, there are no nulls in the magnetic field, and when a = 0, there is a single 
null at 2; = 0. We find that the nulls in the magnetic field are unimportant in this problem 
— rather, zeros of V'^ — A^ are important. The eigenfunctions discussed below (see §6.3p 
show no special behavior at A = 0, but have sharp gradients when \A\ = \V\. In terms of 
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a, \V\ = \A\ at 



z 



±1 



(37) 



a 



(38) 



z 



± 



1 — a 



When a > 0.5, the solutions in eqn. (138!) are no longer in the domain. This means that 
V < A in the entire domain, yielding stability by eqn. f l2^ . Heuristically, when a becomes 
more positive, the strength of the magnetic field in the domain increases until the magnetic 
tension force becomes so strong that all perturbations become stable. 

In the opposite limit, when a becomes very negative, the solutions in eqn. ( ]38|1 
approach z = ±1. For arbitrarily negative a, there is still some region for which V > A. 
Tatsuno & Borland find instability for a as small as -25, and we can prove that there is 
instability for all a < 0.5 when = using the sufficient condition for instability by Chen 
& Morrison described in §11 The explicit computation is messy, but is included in Appendix 



The limit in which a — )• — oo is probably not physically relevant. As the two "almost 
singular" layers approach each other (see eqns. 1371 and [38]) . there are large gradients at the 
boundary of the domain. In this case, the instability probably relies crucially on our choice 
of boundary conditions. Moreover, when stratification is included, the high field strengths 
and large currents corresponding to |a| ^ 1 are destabilizing in themselves, contrary to 
what we assume here. Thus, results in this limit should be viewed as proving a point about 
the Richardson criterion, but are not necessarily physically relevant by themselves. As 
we show in explicit calculations presented below, a does not need to be very negative to 
recover the results described in the infinitely negative case. 
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6.2. Effect of Stratification on Stability 

Our main result is evidence for the following conjecture: There is instability as 
a — 7- — oo, even in the presence of arbitrarily strong density stratification, in violation of the 
Richardson criterion. There does not seem to be any way to prove this claim analytically, as 
there was in the A^^ = case. The sufficient condition for stability presented by Chen and 
Morrison relies crucially on the analytic solution to the eigenvalue equation when fc^ = 0. 
When A^^ 7^ 0, we no longer have an analytic solution to the eigenvalue equation, even when 
/c^ = 0, so there is no extension of the proof. 

Given the assumptions made above, the growth rate c is a function of the following 
parameters: fc^, A^^, and a. We first specialize to the /c^ = case, and then examine the 
more general finite case. 

6.2.1. A;2 = 

For this problem, the unstable area of the (A^^, a) plane is maximized for fc^ = — 
though this is not necessarily true in general (K68). When /c^ = we have c = c(A^^,a). 
We have plotted contours of constant c on the A^^, a plane in Figure [H 

We find instability when A^^ > 1/4, violating the Richardson criterion. It seems that 
given an arbitrarily large value of A^^, there is a sufficiently negative value of a such that 
the fields are unstable. However, as mentioned in §6.2| the extremely negative a case is 
probably strongly affected by the boundary conditions. 

Gravity is stabilizing: the growth rate decreases as A^^ increases. There is stability for 
a < 0.5 by the same arguments as above, and as a becomes more negative, we find larger 
c. Although a stronger magnetic field results in a strong magnetic tension force, and the 
"destabilizing" region in which \V\ > \A\ shrinks for more negative a, we nevertheless find 
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Fig. 1. — Contours of constant c on the a, N'^ plane. The Richardson criterion states that 
the shaded region is stable. The dotted line represents the c = contour — the region below 
this line is unstable. 

stronger instability. We hypothesize that c increases because there is more free energy in 
the magnetic field as a becomes more negative and the magnetic field becomes stronger. As 
a becomes more negative, the instability can tap more free energy from the magnetic field, 
and thus we find a violation of the Richardson criterion. However, note that the stronger 
magnetic field, and corresponding increase in magnetic free energy, is not a sufficient 
condition for instability, as the magnetic field is stable without the presence of shear fiow. 

The contours of constant c are well fit by straight fines. The equation for the boundary 
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between the stable and unstable regimes is 

a = 0.5 - 2.65Ar2. (39) 

Thus, for a < —0.1625, the Richardson criterion is violated. The slopes of the contours 
become steeper as c increases. Although there is instability with arbitrarily large c, this 
does not mean the instability has arbitrarily large growth rate. As mentioned in ^ the 
growth rate is formally zero at k"^ = 0. Thus, to find the growth rate, we need to understand 
the instability at k"^ ^ 0. 



6.2.2. P > 

Although when k^ = Q there is instability for arbitrarily negative a, for every finite 
/c, there is a cutoff for which any a more negative than yields stable profiles due to 
an insurmountable magnetic tension force. Looking at it another way, c always decreases 
as k increases, so for any values a and A^^ which are unstable at k"^ = 0, there is a A; for 
which c = 0. Call this value kcrit{,Oi, N"^). Figure [2] plots kcrit{,Oi, N"^) as a function of a and 
A^^. The point {a,N'^,k) is unstable iff A; < kcrit{o:, N"^)- Although it is possible to find 
instability when fc^ > for profiles which are stable when k"^ = (see §4.31) . this does not 
seem to occur for these classes of profiles. 

Figure [3] plots surfaces of constant u in {a,N'^,k) space. The figure shows that 
a; is a sharply peaked function of k, and that it decreases with increasing A^^. Given 
A^^, ^ 0, there is instability f or only a finite range of a. For A^^ = 0, our results agree 
with iTatsuno fc Dorlandl (120061 ). For sufficiently small k, c is almost constant. Thus, the 
growth rate u = kc is linear in k with slope c. However, as k grows, c begins to decrease. 
There is a maximum growth rate defined by c/logc/rflog/c = — 1, and the growth rate goes 
to zero when c does. The growth rate is 1 - 2 orders of magnitude lower than the typical 



- 26 - 







s -0.5 



1 .0 

1 .5 



max k 



0.100 0.520 0.959 1.559 1.77. 



0.5 




0.05 0.1 0.15 0.2 0.25 



Fig. 2. — The largest k, denoted kcrit, for each a and A^^ which is unstable. The white area 
is stable. 

growth rates of hydrodynamic shear flow instabilities. 

As k"^ increases from zero, the fluid displacement becomes more vertical. Vertical 
perturbations bend field lines, and are subject to a restoring magnetic tension force. Thus, 
it makes sense that the most unstable modes are the horizontal modes characterized by 
k"^ = 0. For some applications, such as stellar interiors (see it is important to consider 
the vertical transport (of angular momentum, etc.) by these modes. In this case, the k"^ = 
mode is irrelevant. One must then consider an optimization problem in which modes with 
too low /c^ have no vertical transport effects, whereas modes with too high /c^ are stable. 



-27- 




Fig. 3. — Surfaces of constant growth rate u in (a, iV^, k) space. The maximum u in this 
range of (a, N"^, k) is given also. 

This argument is only valid assuming that the non-linear evolution is similar over a broad 
range of k"^. A full non-linear simulation for various fc^ is necessary in order to understand 
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the transport properties of these instabihties. 

6.3. Eigenfunctions 

We normahze the eigenfunctions as described in ^ The eigenfunctions all look like the 
example plotted in Figure HI The most salient features are the sharp gradients aX z = ±.47, 
where \V\ = \A\. Notice that the nulls in the magnetic field at a = ±.69 produce no special 
features. 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 



z z 

Fig. 4. — The vertical displacement ^ (left panel) and horizontal displacement = /k 
(right panel), where prime denotes differentiation with respect to 2;, eigenfunctions for a = 
—0.9, A^^ = 0.3 and k = 0.2. The thick solid lines are the real part of the eigenfunctions, 
and the thick dashed lines are the imaginary part of the eigenfunctions. The thin vertical 
dotted lines are at z = ±.47 where \V\ = \A\ and the thin vertical dotdashed lines are at 
z = ±.69, where A = 0. 

7. Constant A with Velocity Profiles Suggested by Kent 



In §4. 31 we summarized Kent's discussion (K68) of velocity profiles which are marginally 
stable in the absence of a magnetic field and destabilized by a small, constant field. In 
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this section we generalize Kent's construction and investigate the stabihty of the resulting 
family of Kent flows. 

The velocity profile is most conveniently specified by the inverse relation z = z{V). 
Note that only invertible velocity profiles, i.e., dV/dz ^ 0, can be specified by this 
inverse relation. When = Q and N'^ = 0, we can use the instability condition by 



Chen fc Morrison 



(Il99ll ) and evaluate the integral in eqn. fl2^ in closed form. This provides 
a transcendental equation for the growth rate. From solving this equation numerically, it 
seems that there exist velocity profiles which are (marginally) stable at = 0, but unstable 
for < y4o < |l^|max- When we increase A^^ from zero, we always find stability when 
A^^ > (maxl^')^/4, but can find instability for all A^^ up to this limit. Our interpretation of 
this result is that the positive energy required to perturb a constant magnetic field triumphs 
over the extra freedom granted by magnetically breaking the frozen- in vorticity constraint. 



7.1. 



First we consider various velocity profiles defined by z = at A;^ = 0. Define 

, , dz 



(40) 



We restrict ourselves to velocity profiles which are marginally stable at A = 0, as they seem 
to be maximally destabilized by magnetic fields. We will first consider velocity profiles with 
walls dX z = ±zo, with the condition that V{-^zq) = ±1. This will simplify the algebra when 
deriving analytic stability results. We will then employ the rescaling symmetry described 
in eqn. ( |T2l) to present numerical results using the normalization zq = 1. 



The condition for marginal stability (jKentlll968l ) is 

u{V)dV _ 
V 



(41) 



30 



where we have assumed that V ranges from —1 to +1 in the domain. Assuming 



z = V + a^V^ + a^V'' + ..., 



(42) 



we have 



10 = 1 + 3a3V^ + 5a5V^ + 



(43) 



so the marginal stabihty condition on the a^-'s is 



E 



J 



1-2 



1. 



(44) 



j>3, odd 



Our construction is a generahzation of K68, who truncated the series in eqn. ( H2|) at 3 
terms. Next we assume there is only one inflexion point, at ^ = 0. This condition implies 
that u cannot have any extrema, so none of the aj are negative. Numerical work suggests 
that the results discussed here hold for velocity profiles with multiple inflexion points, 
so by assuming only one inflexion point, we make the problem much easier, but do not 
qualitatively change the results. 

Now we add a constant magnetic fleld. When /c^ = 0, we have that 



where u}{V) is deflned as in eqn. (140!) . We can rewrite the integral in eqn. (H6|) as 




(45) 



implies instability with growth rate c. If we change variables to V, we flnd 




(46) 




(47) 



The two integrals have equal real parts, so all we need to calculate is 




(48) 
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When specifying uj(y) as a power series in odd powers of V, as in eqn. P3l) . we can evaluate 
the integral by noticing that 

1 c + Aq (c + Ao)3 , , 1 



2 J_^V -c~Ao n-1 n-3 

+ l(c + Ao)Mlog(l-c-Ao)-log(-l-c-Ao)), (49) 

and summing over each term in the power series for uj{V). This gives a transcendental 
condition for stability, instead of the differential condition of eqn. f lTT]) . 

Notice that the location of the walls plays a crucial role in the equation for stability, 
eqn. fl49p . Moving the walls from the Zq where V{zq) = 1 could make the marginally 
stable velocity profiles stable or unstable. Although we will only consider velocity profiles 
which are marginally stable with no magnetic field below, our results do not change 
qualitatively when we add a constant magnetic field to a velocity profile which is stable or 
unstable when Aq = 0. We choose marginally stable velocity profiles because they are more 
clearly destabilized by magnetic fields than unstable velocity profiles, and they are more 
destabilized than stable velocity profiles. 

For the remainder of this paper, we will normalize the problem by setting the walls at 
z = ±1. Under the assumptions that V has only a single inflexion point and is marginally 
stable at A = 0, we numerically find that the most unstable velocity profile ai = and 
N"^ = is given by 

z = V+ 1 + ^ ^ , 50 

\ n J n 

for n odd, when n — t- oo. In this limit, the velocity profile approaches 



V{z) 



+ i, \<z<l 

-\<z<\ (51) 
-i, -l<z<-\ 



For every n odd and greater than three, the velocity in eqn. fl50|) is marginally stable. We 
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plot the velocity profile for n = 5 and n = 41 in Figure [51 Notice that maxV = 1, so the 
Richardson criterion states that N"^ > 1/4 yields stability. 




Fig. 5. — The velocity profile solutions of eqn. (!50|) for n = 5 (solid) and n = 41 (dashed). 

For each n, we can plot c as a function of Aq at = 0. Because we assumed the 
magnetic field is parallel to the velocity, we know there is stability when Aq > Vmax- 
Thus, Vmax sets a natural scale for measuring the magnetic field strength. Figure [6] plots 
c{AQ/Vmax) for n = 'b and n = 41. It seems that as n — )■ oo, the maximum c approaches 
^ 0.125 for Aq ^ 0.65Vmax = 0.325. 

Figure [7] shows an eigenfunction for Aq = 0.65Vmax ~ 0.31, n = 41. Note that it is very 
similar to the eigenfunction for the Tatsuno & Borland profiles in §6.31 
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Fig. 6. — The imaginary part of the eigenvalue c as a function of Ao/Vmax for the velocity 
profiles given by eqn. ( I50l) for n = 5 (solid) and n = 41 (dashed). 

7.2. A^V 

As mentioned in §7.11 the velocity profiles considered here have maxl^' = 1, so the 
Richardson criterion states that A^^ > 1/4 implies stability. As n increases, the maximally 
unstable N"^ increases, but never seems to reach 1/4. Figure [HI shows contours of c as a 
function of A^^ and AQ/Vmax for A; = and n = 41. Although there is instability for A^^ very 
close to 1/4, we find stability at A^^ = 0.25. It seems that the Richardson criterion is not 
violated when adding a constant magnetic field to this class of velocity profiles. 
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-1.0 -0.5 0.0 0.5 1.0 

z 

Fig. 7. — Eigenfunction for A = O.GVmax ~ 0.31 and velocity given by eqn. ( 150|) for n = 41. 
The thick sohd hne is the real part of the eigenfunction, and the thick dashed line is the 
imaginary part of the eigenfunction. The vertical dotted lines denote the points where 
\V\ = \A\, at z = ±0.40. 

Figure [9] shows a typical eigenfunction. As with the velocity and magnetic field profiles 
considered in ^ there are sharp gradients when \V\ = \A\. Unlike the eigenfunctions 
considered above, the real part of this eigenfunction is close to zero at the origin. 

The constant magnetic field case is very different from the parabolic case because 
there is no violation of the Richardson Criterion. We can understand result heuristically by 
noting that a constant magnetic field cannot increase the free energy of the perturbation. 
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Fig. 8. — Contours of c as a function of A^^ and Ao/Vmax for = and the velocity profile 
given by eqn. ( I50l) with n = 41. The white area is stable. 

and thus cannot render a velocity profile with A^^ > maxV"'^/4 unstable. Although there is 
no energy principle in the presence of shear flow, one can show that a sufficient condition 
for stability is that t he energy of a perturbation i s positive, i.e. F(^) ■ ^ > 0, where F(^) 



is the force operator ( iFrieman fc Rotenberg 



1960l ). A constant magnetic field contributes 



+ |QP to the energy of a perturbation, where Q = V x x B). Thus, a constant magnetic 
field always increases the energy of a perturbation. 

However, in §7.H we describe an entire class of velocity profiles which are (marginally) 
stable at y4o = 0, but unstable for Aq > 0. Our interpretation of the destabilized is as 
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Fig. 9. — Eigenfunction for A = O.GVmax ~ 0.31 and velocity given by eqn. ( 150|) for n = 41, 
with A^^ = 0.225, A; = 0. The thick sohd hne is the real part of the eigenfunction, and the 
thick dotted line is the imaginary part of the eigenfunction. The vertical dotdashed lines are 
at ±0.31, where \V\ = \A\. 



1960|), 



follows. An unstable perturbation must have negative energy ( iFrieman fc RotenbergI 
but this is only a necessary condition for instability. Thus, perturbations to the velocity 
profiles considered in §7. II have negative energy, but are still stable. For a sufficiently small 
magnetic field {Aq < Vmax), the increase in energy of the perturbation from the magnetic 
field can be overcome by a negative contribution from the shear fiow, so the total energy of 
the perturbation is negative and there could be instability. 
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Because the Richardson criterion can be understood from energetic arguments (see 
§4.2p . one could assume that when A^^ > \^'^/4 in the entire domain that the energy is 
necessarily positive. Then the addition of a constant magnetic field only further increase the 
energy of the perturbation, preventing instability. This is a rather considerable assumption, 
so this argument is best viewed as a heuristic. 



8. Application to Astrophysical Systems 

We have studied shear flow instability in stably stratified media for flow profiles which 
would be stable in the absence of a magnetic field and shown that Richardson's criterion for 
buoyancy stabilization can be violated, provided that the magnetic field is inhomogeneous. 
In this section we briefly discuss astrophysical applications. 

First, some general considerations. Our analysis holds when the flow and field are 
perpendicular to gravity. We ignored the effect of the magnetic field on the density 
stratification, thereby precluding any instabilities associated with magnetic buoyancy. 
Thus, our work applies primarily to situations in which the field is not too strong and its 
scale height is not much less than the pressure scale height. Thus, although we gave an 
example in ^ of a system that can be unstable at arbitrarily large Ri, instability at large 
Ri required in that case that the flow be sub-Alfvenic in most of the domain and that 
the magnetic scale length be much less than the velocity shear length. In addition to the 
possible introduction of magnetic buoyancy effects, a small magnetic scale height relative 
to the velocity scale height requires that the magnetic Prandtl number Pm — the ratio of 
viscous to magnetic diffusivity — be much greater than unity, opposite to the situation 
in dense plasmas such as stellar interiors. Bearing these things in mind, there is probably 
a practical upper limit on Ri at which magnetic flelds are destabilizing according to the 
mechanism discussed here. 
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It is useful to cast Ri in a form which allows its magnitude to be estimated. We 
introduce a buoyancy parameter fhu in terms of which N"^ can be written in terms of the 
local gravity and pressure scale height as 

= hu^, (52) 

where g and Hp are the local gravity and density scale height, respectively; in the Boussinesq 
approximation, = 1. Specializing to the case that y is a rotational velocity, we introduce 
the velocity scale height H.^ by V' = V/ and a breakup parameter /^r by 

i^'r = /..^, (53) 

V 

where r is the distance from the rotation axis. Using eqns. f l52|) and fl53|) . Ri can be written 

as 

Ri = l^^^. (54) 
thr Hp r 

In stably stratified systems with uniform composition, is generally 0{1), while a 
molecular weight gradient can render ^ 1. Except for systems rotating near breakup, 
fbr ^ 1- Typically, H^ exceeds the geometric width of a shear layer because V changes by 
only a fraction of itself. Thus, although the second and third ratios on the right hand side of 
eqn. are below unity, they are generally not enough to offset fbu/ fbr, and Ri ^ 1. One 
exception to these considerations occurs near the boundaries of convection zones, where A^^ 
crosses through zero. Thus, a thin layer on the stably stratified side of the boundary could 
be magnetically destabilized even if Ri > 1/4. 

The expectation that Ri 3> 1 in the stably stratified portions of stellar interiors is 
borne out by examination of stellar models. First, we consider the Sun. Helioseismology has 
revealed a thin shear layer, known as the tach ocline, below the base of the solar convection 



zone, which is thought to lie at 0.713 Rc:^ (see (Goug. 



at O.7OOi?0 from Gough and V from 



Schatzman et al. 



2007h 



br a review). If we take N'^ 



(|2000[ ). we find that at the equator 
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Ri = 6400 and ~ 10 ^. In other words, even very close to the base of the convection 
zone Ri is quite large, and increases with depth from the value given here. 

We also evaluated Ri in an evolutionary sequence of models of massive, rotating stars 
generously provided to us by G. Meynet. The initial mass is 20 Mq (which decreases due to 
mass loss) and the initial surface rotation period is about 1.2 d. When the star first reaches 
the main sequence, the core is convective and the envelope is radiative. As hydrogen is 
exhausted in the core, strong nonhomologous contraction spins up the core and creates 
strong shear layers, which tends to reduce Ri. At the same time, steep negative molecular 
weight gradients increase fbu- We find that in the bulk of the interior, Ri is between 10^ and 
10^. In the models, the boundaries of convection zones (which form in association with shell 
burning) actually show spikes in Ri. This is because VL is set to a constant in convection 
zones, due to efficient turbulent mixing. Thus, although there is probably a thin layer in 
which Ri drops to small values, it cannot be evaluated from these models. 

These estimates suggest that destabilization of stellar rotation profiles by weak 
magnetic fields is likely to occur only in thin layers outside convection zones. However, 
the tendency for such fields to destabilize a system may be important even when physical 
processes neglected by our analysis are included. Chief among the m is thermal diffusion. 



I974J ) and leads to a larger 



which can suppress the stabilizing effects of buoyancy (iZahnI 
critical Ri to guarantee stabilization. Whether this carries over our analysis is a topic for 
future study. 

The instability could conceivably also operate on poloidal flows. However, because such 
flows are generally slow compared with rotation, their Ri tends to be even larger than Ri 
for rotation. And because rotational shear tends to make the magnetic field predominantly 
toroidal, magnetic effects on the stability of poloidal flow are probably weak. 

Similar considerations hold for accretion disks. The vertical shear in a Keplerian disk 
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of thickness H is smaller than the radial shear by a factor of H/r. If the radial inflow 
velocity is a function of height, its shear could be large, but the magnetic field is expected 
to be predominantly toroidal. Therefore, this instability is probably not critically important 
for either rotation or radial flow in disks. 



9. Conclusion 

Turbulence is a key ingredient in the transport of chemical species, entropy, angular 
momentum, and magnetic flux in astrophysical settings. Shear flows, which are driven 
almost ubiquitously in nature, can become turbulent through instability. 

In this paper we have considered ideal instabilities of magnetized shear flows in stably 
stratified systems. In the absence of magnetic fields, the Richardson criterion provides 
a necessary condition for instability based on comparing the kinetic energy released by 
vertical interchange of fluid elements to the potential energy required to displace them. 
The Richardson criterion is often assumed to set the ideal stability boundary for shear flow 
instabilities in stratified media such as stars and accretion disks. The main result of this 
paper is that the Richardson criterion is no longer valid when inhomogeneous magnetic 
fields are included: because such fields carry free energy, buoyancy forces must be stronger 
to stabilize the system. We have provided an exa mple by adding density stratification to 



the fields described by 



Tatsuno fc Dorlandl (120061 ). These fields can be viewed as a local 



approximation of any shear flow in the presence of a magnetic extremum. The system has 
the interesting property that the flow is neutrally stable in the absence of the magnetic 
field, but unstable in its presence. Solving the eigenvalue problem in eqn. ( ITTl) . we find 
unstable modes for arbitrarily large A^^, provided the magnetic field is sufficiently strong. 
Even for magnetic fields yielding Alfven velocities comparable to fiow velocities, we find 
violation of the Richardson criterion. Thus, when considering the ideal stability of a plasma 
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shear flow in a stratified medium, it is not sufficient to consider the Richardson criterion. 

We were unable to find an example in which a constant magnetic field leads to 
violation of th e Richardso n criterion. We extended and analyzed a class of velocity profiles 



considered by 



Kent 



(119681 ). which were shown to be destabilized by a constant magnetic 
field. Although we were able to destabilize the flows when A^^ = 0, and the fastest growing 
modes have moderately strong magnetic fields, when A^^ > F'^/4, we always found stability. 
We provided two heuristics for understanding the destabilization due to magnetic fields. 
An inhomogeneous magnetic field provides a free energy source which can be tapped by 
an instability. Thus, while a homogeneous magnetic field can be destabilizing because 
vorticity is no longer frozen into the flow, allowing new unstable plasma motions, only 
an inhomogeneous field can provide the source of energy needed to violate Richardson's 
criterion. 

We briefly applied our results to the solar tachocline and to high mass, rapidly rotating 
stars. In the bulk of the tachocline, Ri is very large because the Sun rotates slowly. Very 
near the boundary of the convection zone, Ri drops because A^^ is passing through zero. A 
similar situation holds, for different reason, in high mass stars. Although these stars rotate 
rapidly, the regions of strong shear coincide with regions of strong, stabilizing, molecular 
weight gradient. This keeps Ri large, except near convection zone boundaries. Thus, in 
stars, the destabilization of stratified shear flow by magnetic flelds is most likely to occur in 
thin regions on the stable side of convection zone boundaries. If the weakening of buoyancy 
by thermal diffusion destabilizes magnetized flow in the same way as unmagnetized flow, 
the unstable region could be much larger, however. 

Our 2D slab model is not a realistic geometry for many applications. The introduction 
of additional terms, such as curvature terms from toroidal geometry or the centrifugal 
force for rotation, probably changes our results quantitatively, but not qualitatively. The 
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Boussinesq approximation could also be relaxed to allow more realistic density profiles 
and other physics. Inclusion of diffusive effects would allow us to consider non-ideal 
instabilities, including the secular shear instability. For many applications, the non-linear 
phase and saturation of these instabilities is also important for determining effects such 
as angular momentum transport. These considerations should be investigated further to 
better understand the nature of magneto-shear instabilities in a stratified medium. 
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A. Instability of V = z, A = {1 — a)z'^ + a when a < 0.5 

We will prove that the velocity and magnetic field profiles considered in ^ ^ = z, 
A= {1 — ajz^ + a, are unstable when a < 0.5. In §4.31 we described the following sufficient 
condition for instability at /c^ = by Chen and Morrison (eqn. flMl) ): If 

as e — )■ 0, then there is instability. We can factor the denominator to get 

1 f^^ dz 1 f^^ dz 

(A2) 



2 Aiy-ie-A) 2 7_i Aiy-ie + A) 
Let us examine how these two integrals are related. Define u = —z. Then 
1 /•+^ dz 1 f-^ du 



2 7_i A{z){V{z) -te + A{z)) 2 J^^ A{z){V{z) - ie + A{z)) 



(A3) 
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_i du 

~2y_i A{u){-V{u)-i€ + A{u)) 
1 /•+^ rfw 



2y_i A{u){y{u) + ie- A{u)y 

which has the same real part as the first integral, but opposite imaginary part. Thus, we 
need only check that 

as e — 7> to prove instability. Integrals of this form can be evaluated in closed form, but 
must first be factored. To simplify the algebra, we reduce the degree of the polynomial in 
the denominator through partial fractions. 



3? 



r dz dz dz 

J_, A{V-te-A)^^J_, A{V-ie)^^J_, {y - ie){V - le - A) ^^^^ 



The first integral gives no contribution because multiplying hy V -\-ie in the numerator and 
denominator shows that the real part is odd and integrates to zero. Thus, we need only 
evaluate the second integral. 

We can integrate the remaining part by brute force, i.e. using Mathematica. Assuming 

e > 0, Mathematica gives 

/ iy-^e){V-^e-A) = ' Aia-l + ae^) V^' (D ^^'^ 

- log (e^ + (-1 + zfz" - 2a{-l + zfz{l + + + z^f) 

Ail - 2i{l - a)e) ( -1 + 2(1- a)z 

+ ^ ^ ' ' arctan ' ^ ' 



•y/-l - 4q;2 + Aie - Aia{i + e) y -yZ-l - Aa'^ + Aie - 4ia{i + e) 

f (—1 + z)(z — az — a)\ , /o ox 
+ 2i arctan I ^ ^ J + 2 log {e^ + z^) . 

Notice that the prefactor has the opposite sign as a. The term on the first line is imaginary, 
so we do not need to consider it. In the logarithm on the second line, the third and fourth 
terms which are at 2; = ±1. On the last line, the first term is imaginary and the second 
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term is even, so neither contribute to the integral. Thus, if 
1 



4(q; - e2 + ae2) 



log ( ^ ) m 




4(l-2z(l-a)e) 
+ 3f? , ^ \ ^ ^ arctan 

- 4a2 + Aie - Aia{i + e) 

^ 4(1-22(1 -a)e) . . , ^ 

- ^ . arctan ^ 1 > 

- 4a2 ^ _ 4^Q,(^^ _^ y _ 4(^2 ^ 4^g _ 4^^(5; _^ 

for a particular a as e — )• 0, then the profiles for that a are unstable. The e for which the 
RHS of eqn. f lA7p equals zero is the growth rate of the instability. Thus, this relation gives 
a transcendental equation for the growth rate, which is significantly easier to solve than the 
differential eigenvalue problem given in ^ 

As e — 7- 0, the logarithm term diverges and is positive. However, when z = +1, the 
arctan term also diverges, approaching —ioo, meaning that the entire term gives a negative 
divergent contribution. We need to see which diverges faster. The argument of the z = +1 
arctan term is 

1 - 2a . 1 - 2a 



^-1 - 4a2 + Aie - Aia{i + e) ^/4a^ - 4a + 1 - 42e(l - a) 



-I 1 



4a2 - 4a + 1 
14ie(l - 
2 (l-2a) 



+ (A8) 



In general, arctan(z) is given by 



arctan(z) = i- (log(l — iz) — log(l + iz)) . (A9) 

The divergent part for us is the first term, so 

/ -1-2(1 -a) \ 1 / 2ie(l-a)\ , 
arctan ^ ^ ^ t- log -— ^ . AlO 

If we neglect the e terms which are not in the divergence, we find that the coefficient of 
the log(e) term is —2/(1 — 2a). Thus, only considering the terms in eqn. ( ]A7I) which are 
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divergent as e — > 0, and taking e = except for in the divergence, we are left with 



When q; < 0, we have that —1/Aa > 0, and the first log(e) term dominates, so the whole 
quantity is positive. Thus, we have proven that there is instability for a < 0. When 
0.5 > « > 0, we have — l/4a < 0, but the second logarithm term dominates and is negative, 
again yielding instability. However, when a > 0.5, both divergent terms become positive, 
but —l/Aa < 0, so the quantity is negative as e — > 0, and the profiles are stable. In order 
to show instability at a = 0, we would need to retain more terms in our perturbative 
expansion in e. 




(All) 
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